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SUMMARY 

A description is given of a measurement method, applicable to a polydis- 
persion of particles, in which the intensity of scattered light at any angle is 
weighted by a factor proportional to that angle. Determination is then made of 
four angles at which the weighted intensity is four fractions of the maximum 
intensity. These yield four characteristic diameters: the diameters of the 
.volume/area mean (D32 the Sauter mean) and the volume/diameter mean (D31); the 
diameters at cumulative volume fractions of 0.5 (Dyo.5 the volume median) and 
0.75 (D v q 75)- They also yield the volume dispersion of diameters. Mie scat- 
ter computations show that an average diameter less than three micrometers 
cannot be accurately measured. The results are relatively insensitive to 
extraneous background light and to the nature of the diameter distribution. 

Also described is an experimental method of verifying the conclusions by using 
two microscope slides coated with polystyrene microspheres to simulate the par- 
ticles and the background. 


INTRODUCTION 

This paper deals with the problem of determining by optical methods (a) 
the mean diameter of a large aggregate of particles such as occurs in fuel 
atomization studies or in aircraft icing studies and (b) a measure of the dis- 
persion of particle diameters. The particles are produced by atomization of a 
continuously flowing liquid. The number of particles in the field of view is 
usually greater than 100. 

When the particle diameter is greater than the wavelength of radiation, 
the diameter can be inferred from the angular distribution of light scattered 
in the forward direction from a collimated beam. 

Reference 1 shows that a useful measure of particle diameter is derivable 
from the shape of the center lobe of the intensity distribution of forward- 
scattered radiation. The intensity may be calculated by the Fraunhofer dif- 
fraction formula instead of the more exact Mie formula. Although there is a 
substantial discrepancy between the two formulas in the calculation of absolute 
intensity when the particle diameter approaches the wavelength of the radia- 
tion, the absolute value of intensity can be eliminated by measuring the ratio 
of intensities at two fixed scattering angles. The lower limit of measurable 
particle diameter then becomes of the order of the wavelength of the radiation. 
The upper limit of measurable particle diameter is set by the increasing diffi- 
culty of measuring intensities as the angular width of the center lobe becomes 
smaller for a larger particle diameter. When measuring particles larger than 
10 pm in diameter, the intensity in the center lobe is due principally to 
Fraunhofer diffraction rather than to the additional processes of reflection, 
refraction, and absorption that are included in Mie's analysis. 



Reference 2 , an application of reference 1 , measured the angles, one fixed 
and one variable, where there exists a fixed ratio of two intensities. The 
methods of references 1 and 2 were compared by reference 3 . However, both 
methods are limited in accuracy by extraneous background light. The method 
presented here enables the two angles of measurement to be chosen in such man- 
ner as to reduce the effect of extraneous light. This method is similar to 
that used by references 4 and 5 to determine particle diameter distribution, 
except that in those references, the effect of extraneous light was not consid- 
ered; more recently, the effect of extraneous light was treated by references 3 
and 6, and it is treated in detail in this report. 

When there is a polydispersion of particles, it is generally recognized 
that the Sauter mean diameter ( D32 ) is a good parameter to use in studying fuel 
combustion rates, and the volume median diameter (D v q 5) is a good parameter to 
use in studying icing conditions for aircraft in clouds and rain. The polydis- 
persion may be characterized by its dispersion (the span of particle diameters) 
and by its shape as given by the analytic form of the distribution function. 
References 7 and 8 show that the light-scatter curve determines D32 nearly 
independently of the particle diameter distribution function. In the present 
work, this independence is confirmed for a narrow range of dispersions, but, 
contrary to references 7 and 8, there is found to be a strong dependence on the 
dispersion of particle sizes. To minimize this dependence, we find a point on 
the scattered-light intensity curve at which the intensity is a specified frac- 
tion of the maximum intensity. Four intensity ratios are chosen so that the 
diameters D32, D31 , D v q. 5 and D v q. 75 are determined, and these values are 
substantially independent of both the distribution function and its dispersion. 
Furthermore, the dispersion is also determined by these measurements. The 
ratio D v o.75/D v q. 5 may be useful as a measure of the relative amount of large 
diameter particles. This ratio has served to indicate the transition through a 
bimodal volume distribution as air assist in a spray nozzle was increased from 
zero. The ratio D31/D32 may be useful to better characterize the spray when 
there is evaporation. 

This paper also describes a means of simulating the actual physical situa- 
tion by an optical-laboratory experiment. Commercially available microspheres 
are deposited on glass plates to simulate both a spray of droplets and the 
extraneous background light that interferes with the measurement. The methods 
of properly correcting for this background are presented and experimentally 
conf i rmed . 


TERMINOLOGY 

Symbols used in this paper are presented in Appendix A. A further defini- 
tion of some terms will be given here. 

Ordinarily, a large aggregate of liquid droplets will consist of droplets 
of varying diameter D, ranging from a minimum value D& to a maximum value 
D u . In any given volume V 0 , the number of particles whose diameters lie in a 
certain narrow range of diameters is conveniently represented by a number dis- 
tribution function 


dn/dD 
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where dn denotes the number fraction of particles lying in the diameter range 
dD and 


(dn/dD >dD = 1 


Similarly, one may define a volume distri jution function 

dV/dti 

where dV represents the total volume of all particles lying in the diameter 
range dD and 

(dV/dD) ID = V Q 


For reasons of mathematical convenience or of convenience in physical 
analysis, the lower limit may, on occasion, be taken either as zero or as a 
small multiple of the wavelength of light being used. For the same reasons, 
the upper limit D u may be taken as infinity or as some arbitrary finite 
value. 

The aggregate of particles may often be characterized by (a) a single cen- 
tral value of diameter and (b) a single parameter that represents the disper- 
sion of diameters about that central value. In fuel combustion studies, the 
most useful central value has been found to be the Sauter mean diameter D 32 
given by 

u (dn/dD) • D 2 dD 
ft 

In aircraft icing research, the most useful central value has been found to be 
the volume median diameter D v o.s given by 

D 

V JZ » I V °‘ b (dV/dD)dD 
O J n 


The terminology used here follows that of ASTM Standard Practice E799-81 , 
reference 9, which may be consulted for more detail. 


THE SCANNING-SLIT METHOD 
Mechanical Arrangement 

An optical system described in references 7 and 10 for measuring the radi- 
ant flux at any small angle of forward scattering Is shown in figure 1. Light 
from a point source is collimated by a lens, passes through the test section, 
and is then focused to a point in the focal plane by an objective lens. Light 
scattered anywhere in the test section at an angle © intersects the focal 
plane at a distance y from the focal point. Thus, 
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e = y/f 


(1) 


for y/f << 1, where f is the lens focal length. At the focal plane, a plate 
with a pinhole area SA subtends a cone with apex at the lens and with its 
axis in direction e. The radiant power received by the area SA determines 
the irradiance as H = P/SA, power/area. The intensity in direction 9 is 
I = Pf2/sA, power/steradi an . Thus, the intensity and irradiance are related by 
only a constant as 


H = l/f2 (2) 

This conversion is useful when calculating the power available to a detector at 
a given intensity. The scanning pinhole in figure 1 thus transmits the irra- 
diance distribution which is proportional to the intensity distribution at the 
test section by equation (2). 

To scan over a range of the angle 9, one or more pinholes may be part of 
a rotating disc to produce a periodic electrical waveform. Or, an image-type 
detector may be used to obtain nonmechanical scanning. 

An alternative scanning slit assembly that has an area SA proportional 
to angle 9 is shown in figure 2. It consists of a fixed V-shaped aperture 
and a slit that is translated in a radial direction. This produces a weighting 
factor proportional to 9 that multiplies the intensity distribution provided 
by the scanning pinhole alone. (Since the radial translation of the scanning 
slit is mechanically awkward, the slit is usually cut into a motor-driven 
rotating disc; a weighting factor proportional to 9 can be maintained with 
less than one percent error by using the curved aperture shown in fig. 3.) 

The intensity distribution of the scanning pinhole and the scanning-slit 
assembly will be compared. 


Intensity Distributions 

Useful approximations for the intensity distribution of forward-scattered 
light in the center lobe (the lobe with maximum intensity at 9=0) are plot- 
ted in figure 4. The ordinate is the intensity at a scattering angle 9 
divided by the intensity at 9=0, the direction of the incident light. The 
abscissa is a beam-spread parameter (a9) 2 where the particle size parameter is 

a = trD/X. (3) 


with particle diameter D and wavelength X. 

The curve for a monodispersion is given by the Fraunhofer diffraction for- 
mula (ref. 1) as 


1(9) 

1(0) 


2J^ (a sin 9) 


-.2 


a sin 9 


( 4 ) 


where J] is the first-order Bessel function. For 9 < 20°, sin 9 may be 
replaced by 9 in radians. 


4 



The two curves for a polydispersion of particles are based on the Sauter 
mean diameter D32 - This diameter is the ratio of the third moment to the sec- 
ond moment of the particle-number di stri b Jtion . The curves were derived in 
references 7 and 8 using equation (4) with a variety of particle-number distri- 
butions in which no particles exist with 1 diameter larger than a given multi- 
ple of that diameter at which the derivative dn/dD is a maximum. 

The Gaussian curve is a good approximation to the other curves, and it 
will be used for error analysis, because of mathematical convenience. 

Figure 5 shows the Gaussian relative intensity curve of figure 4 plotted 
with abscissa a© instead of (a©) 2 , and ilso another curve in which the ordi- 
nates have first been multiplied by a weighting factor proportional to © (as 
would be produced by the scanning-slit assembly) and then normalized to a maxi- 
mum value of unity for the relative intensity. One significant difference is 
apparent. With the scanning slit the maximum of the curve is moved away from 
a© = 0 to a© - 1.25. If particle diametar is to be inferred from the angle at 
which the weighted intensity is a fixed fraction of its maximum value, only one 
angle must be measured accurately. For e<ample, if the intensity ratio is cho- 
sen to be 0.4, then a© must be 2.6. The angle at the broad maximum is not 
explicitly measured; instead, only the amplitude of the weighted intensity max- 
imum is needed. In this sense, only the one angle at a fraction of the maxi- 
mum is measured. 

Measurement of the maximum and the angle © at a fraction of maximum can 
be made in real time with an electrical waveform peak detector as used in ref- 
erence 2; or, with a small time delay, by computer using one of several commer- 
cially available instruments designed for electrical waveform acquisition, 
storage and processing. The latter method is usually preferred. 

The extraneous light produced by dif c raction at the test section aperture 
is substantially ignored by the type of measurement described above, as shown 
by the following analysis. 


Diffraction by Test- Section Aperture 

Our tests of a forward-scatter optical system have shown that a fundamen- 
tal limitation to measurement accuracy is extraneous light produced by diffrac- 
tion of the light at the test section ape'ture even when the instrument optics 
produce negligible extraneous light by reflection or scatter. 

The intensity distribution described by the Fraunhofer diffraction equa- 
tion (4) applies both to the spray particles and to the test-section aperture. 
The existence of diffraction by the aperture is often ignored or overlooked. 

The methods used in references 11 and 12 eliminate the effect of diffraction 
by the aperture as well as the effect of other sources of interfering light. 
These methods may be considered when extraneous light is a major problem not 
controlled by the present method. 

Because the test-section aperture diameter is much larger than a particle 
diameter, the aperture produces more scattered flux than does a small number of 
particles. This flux has a strong effect on the resultant intensity distribu- 
tion as measured. 
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The Fraunhofer equation gives the intensity as: 


1 ( 9 ) 


hA 4 

16ir 2 


2J,(a6) 


-.2 


a0 


( 5 ) 


where H is the irradiance of a plane wave of radiation at the aperture. 

The Bessel function J] has peak values, the first one beginning at 
a0 = 1.84, that are closely approximated by an asymptotic formula (ref. 13): 

J 1 ( a© ) = (2/ira0) 1/2 COS(a0 - 3ir/4) (6) 


Equation (6), with the cosine factor replaced by unity (its maximum possible 
value), is substituted in equation (5) to yield an approximate envelope of the 
intensity peaks for a© > 1.84. Along this envelope, for a© > 1.84, 


I (©) 



(7) 


Figure 6 shows in solid lines the relative intensity I(©,D)/I(0,do) as a 
function of angle 0, for a uniformly illuminated circular aperture or spheri- 
cal particle of diameter D(when do = 10 000 pm = 1 cm) as measured by a scan- 
ning pinhole. For a© > 1.84, the intensity has been computed by equation (7). 
For a© < 1.84, the intensity has been computed by equation (5). For 
a© < 1.84, the two formulas yield the same numerical value within 2 percent. 

To determine particle size solely from two intensity measurements, one 
makes measurements in that part of the center scattering lobe which corresponds 
to the knee of the curves. For a single particle and a circular aperture, fig- 
ure 6 shows that at the knee of the curves for smaller diameter particles the 
intensity due to the uniform beam is greater than that due to the particle by a 
factor equal to the ratio of the aperture to particle diameter. This is unde- 
sirable because the uniform beam intensity overpowers the particle intensity. 
However, the situation can be improved by reducing the intensity diffracted at 
the aperture by use of aperture apodization, treated in reference 14. 

Alternatively, instead of using a circular aperture with uniform irra- 
diance, the irradiance may be made to decrease with distance from the optical 
axis. A laser beam, with the usual Gaussian irradiance distribution, has this 
property. The laser beam is optically expanded to a desired diameter at the 
test section. The Gaussian beam may be nominally sized by the "diameter" at 
which its irradiance is some fraction of the axial irradiance. In figure 6 the 
dotted lines apply to a Gaussian beam for which the irradiance amplitude ratio 
H/Ho is 0.5 at a diameter of 1.0 cm. For this beam, an aperture of diameter 
d = 1, 2, or 2.58 cm would produce H/Ho equal to 0.5, 0.0625, or 0.01 respec- 
tively. The improvement shown by the dotted line curves for these three aper- 
ture diameters is substantial and rapid at small scatter angles, after which 
an asymptotic slope l/©3 is followed. The computation was done with the 
method of reference 15; that is, by a summation of Bessel functions of integer 
order instead of by analytic integration. 

The gains obtained by apodization are also possible by spatial filtering 
(references 3, 16, and 17) but with a more complicated optical system. 
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For an unweighted measure of the spra/, with uniform irradiance, and to 
minimize the weighting by a Gaussian beam, the beam in the test section should 
be rectangular. Thus, a circular aperture may be reduced by a wide slit with 
the long dimension oriented perpendicular to the direction of the spray veloci- 
ty. The slit edge diffraction will be small in the direction nearly parallel 
to the slit edge. This modification is recommended. 

Figure 6 shows that the intensity ratio would be very low if one sought to 
measure a single small-diameter particle. However, in the application consid- 
ered here, there is a constant volume flow rate of liquid, and the number of 
small-diameter particles becomes very large. 

The transmission factor is related tc the number of particles, and there 
is a gain in intensity of the scattered light when. N particles are present. 

If N particles of diameter Dp are distributed in a test section of diameter 
. Db, the transmission factor (ref. 10) is: 

( -KN-rrD^/4\ 

x(N) = exp — (8) 

V J 

where K is the scattering coefficient of a particle, equal to 2 for a parti- 
cle diameter larger than the wavelength. When K = 2, some of the flux inter- 
cepted by the cross-sectional area is scattered by reflection and refraction 
through the particle, and an equal amount of flux is Fraunhofer diffracted on 
passing outside the particle. Thus, the diffracted radiation from N parti- 
cles is equal to one-half the flux removed from the beam, and is a fraction 
[1 - t(N)]/ 2 of the incident light. The gain over the light intensity from 
one particle is thus: 


1 - x(N) 1 
1 - x( 1 ) = 2 



(9) 


where, since x(l) is near unity, the first two terms of the series approxima- 
tion of x are adequate to replace the denominator of equation (9), and K 
has been set equal to 2. 

The intensity produced by the test section aperture is reduced by the 
transmission factor x(N). Thus, the overall effective gain in intensity is 
the gain in equation (9) times l/x(N): 


_ [1 - x(N ] 

2x(N) ^D p J 


(10) 


If a single droplet were replaced by N droplets having the same total 
volume, the number N of droplets would vary inversely with the cube of drop- 
let diameter. For example, if one drople 1 : were 10 3 pm in diameter there would 
be 10 3 or 10^ droplets of diameter 100 or 10 pm, respectively. Table I shows 
how the values of x and G would changv with diameter. The intensity ratio 
for this example is shown by the dashed 1 nes in figure 6. The favorable 
results are: 
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(1) A large reduction in the range of variation of relative intensity with 
droplet diameter. 

(2) At the knee of the curve, the effect of the aperture, for uniform 
illumination, is negligible for droplet diameters less than 10^ ^m. For Gaus- 
sian illumination with H/Ho = 0.0625, the effect of the aperture is small for 
droplet diameters less the 10^ pm (0.1 the beam diameter). Thus, a large range 
of droplet diameters can be measured. 

An unfavorable result is that the transmission factor in table I is less 
than 0.2 for 10-pm diameter particles. Multiple scatter at small values of x 
can cause a smaller measured mean diameter as shown by reference 18. 

To apply to the scanning-slit method, the ordinate of the solid-line 
curves of figure 6 is multiplied by a weighting factor of 1000 O to produce 
the solid-line curves of figure 7. A result evident from figure 7 is that the 
measured intensity due to particles changes very little with particle diameter. 
This helps measurement over a large range of diameters. 


Effects of the Form of the Distribution Function and of the Dispersion 

The intensity distribution with a scanning pinhole (fig. 4), was shown by 
references 7 and 8 to be relatively insensitive to the form of the particle 
number distribution function. It is necessary to find if this is also true for 
the intensity distribution with the scanning slit (fig. 5). In the following 
discussion particle volume distribution will be used instead of number distri- 
bution because it is a more widely accepted measure of spray produced by noz- 
zles. Four forms of distribution functions as described by reference 19 will 
be used: Rosi n-Rammler (RR), Nukiyama-Tanasawa (NT), log normal (LN), and 

Upper Limit Distribution Function (ULDF). The ULDF will be used with two val- 
ues of the upper limit. Two more functions not so well known will also be 
used: Wei bull, reference 20 and Rinkes, reference 21. 

A cross-section of the spray near a nozzle has a mean particle diameter on 
axis, that may progress to a substantially different mean diameter off axis. 
Such a spray has a volume distribution different from the six functions to be 
used. This distribution will not be considered because it is more correctly 
determined by measurements with a narrow laser beam at a number of distances 
from the spray axis. These measurements may then be inverted to give the mean 
diameter as a function of distance from the spray axis. An appropriate inver- 
sion method is described by reference 22. 

The spread of any one distribution curve will be characterized by the dis- 
persion A v , where 


A v = 2 (D v q.9 - D v q. l )/(D v o.g + D v q ] ) (11) 

and where D v o.i(D v o.g) = particle diameter such that one tenth (nine tenths) 
of the total volume of all drops is composed of drops smaller in diameter than 
DvO.l^vO^- 

If D represents the diameter at which the distribution is a maximum and 
if the distribution functions are normalized so that this maximum is unity, the 
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six distribution functions appear as plotled in figure 8, for an arbitrary 
diameter D = 50 pm, and a dispersion A v = 1.35. This dispersion represents 
the widest usually produced by a spray, ihese distributions are considered to 
bound the actual distribution produced by a spray nozzle. The curves are 
represented by the following formulas in vhich 6 = distribution parameter, and 
x=D/D; z=(dV/dD)/(dV/dD)o 

(1) Rosin-Rammler, RR, 


z = e A-B ( 6-l 


( 12 . 1 ) 


where 

A = (6-D/6 
B = Ax s 

(2) Nukiyama-Tanasawa, NT, 


z 


e A -*x5 


( 12 . 2 ) 


where 
A = 5/6 
B = Ax s 

The RR and NT functions are identical whei 6-1=5. 

(3) Log Normal, LN, 


z 


X 


2 2 
-1 (/ -B ) 


where 


A = 0.56 B * 6 • 1 n x - A (12.3) 

(4) Upper Limit Distribution Function, ULDF, where D u is the largest 
particle diameter, 


z 


2 

(A -C ) 
e 


/(Bx) 


(12.4) 


where 


x^ = Du/D 



B = (x u - x)/(x u _ 1 ) 


C = 6»ln(Bxe ^) 
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(5) Heibul 1 , 


A-B 6+2 
e x 


(12.5) 


2 = 


where 

A = (6 + 2)/S 
B = Ax 5 

The Weibull is identical to the RR when 6+2 is replaced by 6-1. The 
Weibull is identical to the NT when 6+2=5. 

( 6 ) Rinkes, 


z = (A/x 2 ) e - B(1 - x > /x 


( 12 . 6 ) 


where 

A = (6 - x ) / ( 6 - 1) 

B = (56 - 4) / ( 6 - 1) 

6 - D u / D = x u 

The quantity D u is the upper limit of particle diameter. The parameter 6 
in these distribution functions governs the spread in particle diameters. The 
relation between the dispersion A v and the 6 's in the various functions 
given by equations 12.1 to 12.6 is shown in figure 9. 

As shown by reference 19 the RR function has unreasonable values of mean 
diameter. The cause for this is shown by the relative number distribution 
plotted in figure 10. The curves have been normalized so that they have a max- 
imum of unity. The RR and Weibull functions go to infinity at 2 ero diameter 
for this large dispersion; consequently they have been truncated, so that they 
remain at a value of 2 ero for all diameters D less than D& = 0.1 pm 
(D ft /D = 0.005). 

The effect of the various distributions on various definitions of mean 
diameter Dpq is shown in figure 11 for all combinations of p and q for 
which p + q\ 5. Three truncation points have been assumed for the RR func- 
tion; only the truncation at Dj^ = 2 pm (D^/D = 0.025) appears to yield realis- 
tic values of mean diameter, comparable to those yielded by the other 
di stribution_functions . The RR function will therefore be used with trunca- 
tion at Dq/D = 0.025. The Wej_bull function will be accepted, with truncation 
as shown in figure 10, at Dq,/D = 0.002. 

Figure 12(a) shows the scattered-light-intensity distributions for the 
seven distributions assuming a constant value of D v q. 5 , and three values of the 
dispersion A v . For any one value of dispersion, all Intensity distributions are 
sufficiently alike so that only the envelope of each band of curves is shown. 

For Ay = 0.1, the band is so narrow that only a single curve need be drawn. 
Distributions are shown both for a scanning pinhole (with ordinate I/I max ) anc * 
for a scanning slit (with ordinate I w /Iw,max>- Figure 12 (b) shows a similar 
set of distributions assuming a constant value of D 32 , figure 12 (c) is for a 
constant value of Dyo.75, and figure 12 (d) is for a constant value of D 3 ] . 
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The abscissa in figure 12 is the beam spread parameter <x6. The maximum 
spread in the value of the abscissa at any ordinate, after all distribution 
functions and all dispersions are considered, may be represented by the maxi- 
mum fractional deviation A(a0)/(a0). 

In figure 12, the value of the ordinate at which there is smallest depend- 
ence on the dispersion and on the distribution function is determined as fol- 
lows: A median curve is drawn for each pair of envelopes (for A v = 0.1, the 

curve drawn is already the median); the point sought is taken as the center of 
area of the triangle formed by the three points of intersection of the three 
median curves. (For greater accuracy in determining the centroid, the curves 
were redrawn to greatly enlarged scale in the vicinity of the intersections; 
however, these enlarged curves are not sho^n here.) In figure 12(c), for the 
scanning slit, the median curves are closest together on the left of the maxi- 
mum. The following conclusions result. 

(1) The coordinates of the points of the smallest dependence are as listed 
in table II, which also lists the maximum spread A(a0)/(a0) = A0/0 in the 
abscissas of those points. 

(2) If intensity measurements are made at these points, the value of 

D v o 5 and D v q 75 will be less dependent on the choice of volume distribution 
function than will the value of D32. This conclusion holds regardless of 
whether the scanning slit or the scanning pinhole is used. 

(3) The use of the scanning slit requires measurements of A0 at a much 
larger fraction of the maximum intensity than does the use of the scanning pin- 
hole. This improves accuracy of measurement. 

(4) Away from the selected values of A(a0) listed in (1), there is a 
strong dependence on particle volume dispersion. This result does not agree 
with the conclusions of references 7 and 8, because the dispersion was not var- 
ied over a wide range, in those reports. 

(5) The envelopes of figure 12 are bounded principally by the RR on the 
lower part and by the Rinkes on the upper part, following the progression of 
D oa in figure 11. When the distribution function is not known, one may choose 
a typi cal function near to the center of the group like the NT or the ULDF with 
D u /D = 10. 

The abscissas of all measurement points listed under 1(a) to 1(g) lie 
between 0.9 and 3.3. In order to clarify the dependence on dispersion, the 
curves of figures 12(a) to (d) have been redrawn in figure 12(e) so that they 
pass through a common intersection point at an abscissa of a 0 = 2.5 with an 
intensity ratio of 0.4 for the scanning slit and of 0.04 for the scanning pin- 
hole. The common intersection point was produced by letting the diameter D 
of each function deviate from 50 pm. The envelopes are thereby condensed to 
show explicitly the effect of dispersion and of the function (principally skew- 
ness of the function). At the intersection, the slope of the envelopes shows a 
significant change with dispersion. The width A(a0> of an envelope even at 
the largest dispersion is narrow relative to the change of a0 with disper- 
sion. Because the envelope contains the curves of all functions, this method 
is insensitive to the distribution function. 
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Dispersion of Particle Volume Distributions 


For any one curve, we shall define, for the scanning slit, a slope parame- 
ter (sf) as 


( I/X max) 


32 


s e = 


r-Xo, 


-0.42 


<a9) 


32 


(« 0 ) 


(a0) 


- 1 


32 


v0.5 


(a©) 


- 1 


v0.5 


(13) 


where 

<W max^v 0.5 = O - 73 = optimum intensity ratio for measurement of D v q .5 
(Iw/Iw, max^32 = O - 42 = optimum intensity ratio for measurement of D 32 

Qv 0.5 = an 9 le at wh1ch < V^.raaxJvO.S = °- 73 
e 32 = angle at which < I w /I w ,max>32 = °- 42 


Using the defintion of A v given by equation (11), the plot in figure 13 
shows the mean value of s(A v ) and, by the bars at the various points, the 
spread in the values of s for the seven distribution function shown in fig- 
ure 10. Figure 13 shows an excellent correlation between Sf and A v that 
is substantially independent of the volume distribution function. Thus, fig- 
ure 13 may be used to determine the volume dispersion A v . Figure 9 may then 
be used to determine the parameter & for any of the functions in equa- 
tions 12.1 to 12.6 if these are of interest. The average deviation (a.d.) of 
Individual points from the curve drawn in figure 13 is 0.4 times what the a.d. 
would be if 2 D v q .5 were used in the denominator of equation ( 11 ) instead of 

(DvO. 9 + D v o.l>- 


Procedure for Treatment of Experimental Data 


The preceding analysis leads to the following analytical treatment of 
experimental data obtained by the scanning-slit apparatus. (Note that this 
apparatus yields I w as a directly measured quantity.) 

(1) By varying angle ©, the maximum weighted Intensity I wmax is found 
and recorded. 

(2) By varying angle ©, the angles 0 v q 5 and ©32 are found at which 
I w has a value given by I W) max x (ordinate listed in table II). These angles 
are recorded. 


(3) The slope factor s is computed by equation (13) and the dispersion 
Ay then determined from figure 13. 


(4) The diameters D v q .5 and D 32 are computed from the measured angles 
©v0.5 anc * e 32 ancl t he abscissa values listed in table II. 
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(5) Operations like steps 2 and 4 may also be used to compute D v q .75 
and D 31 , if they are of interest. 

An alternative procedure that leads drectly to the relation between par- 
ticle diameter and angle, is possible if there are available several slides, 
each of which has a monodispersion of representative particles in the diameter 
range of interest. With apparatus like that to be described later under 
Experimental Simulations, such calibration can be more accurate than step 4, 
above . 

A correction of the measured slope fo ' possible instrument error may 
improve accuracy. When diameter calibration is made with microspheres having 
dispersion near 0.16, the measured slope may be compared with Sf in figure 13 
and equation (13) to determine a correction. 

Two additional refinements are possib e. One is a correction for the case 
where particle diameters are so small that Mie's treatment of diffraction phe- 
nomena is significantly superior to Fraunhofer's method. The other is a cor- 
rection for multiple dispersion when parti : 1 e density is high. These two 
corrections will now be considered. 


Corrected Computation Based on Mie Diffraction 

Up to this point, the dispersion and 'our mean diameters have been deter- 
mined analytically assuming Fraunhofer diffraction. However, for small diame- 
ter, Mie's treatment of diffraction is necessary to obtain a true intensity 
distribution. Such a computation has been made for the scanning-slit method, 
and the results expressed in the form of corrections to be applied (itera- 
tively, if desired) 

( 1 ) To the value of the ordinate to be entered into figure 13, in order 
to find a closer value of A v , and 

(2) To the value of the originally computed mean diameter, in order to 
find a closer value. 

The procedure of the computation followed that which was used to obtain 
figures 12 and 13. For the seven distributions previously assumed, and for 
five values of dispersion, the parameters (a 0 ) at the four ordinates listed in 
table II, and the corresponding slope factor s computed from equation (13), 
were determined. Calling these quantities (a©)^ and sy, respectively, and 
calling the respective original values (aG)p and sp, figures 14 and 15 result 
for Mie scatter with water. 

The corrected procedure then becomes: 

(1) Use the original procedure (steps 1 through 4) to obtain first esti- 
mates of dispersion A v and mean diameters D v q .5 and D 32 . 

(2) Use the measured slope sm and the equation in figure 14 to find 
sp/SM- Then compute an equivalent s as 

s = (Sp/SM x s M 

and enter this value of s into figure 12 to find a corrected A v . 
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(3) Use figure 15 to find a corrected mean diameter from a corrected value 
of 9 given by 


0 = (G^/Qp from fig. 15) x (measured 9) 

This last step is not necessary if an experimental calibration is available 
with microspheres in the range 3 to 12 pm. 

The graphs of figures 14 and 15 represent computed points joined by 
straight lines. Each point is the median among the seven assumed distribu- 
tions. Vertical bars at the points represent the spread among the seven dis- 
tributions. At the right-hand end of figure 14 are shown the original values 
of sp and their associated spreads when Fraunhofer diffraction had been 
assumed. Similarly, the spreads among the seven distributions are shown in 
figure 15, along with the value of A0p/9p from table II. 

The following deductions are noteworthy: 

(1) Corrections become significant only when mean diameter is less than 

12 pm. 

(2) For A v > 0.5 and D v q .5 < 3 pm, the differences among the curves are 
comparable to or smaller than the spread for any one curve. Consequently, the 
lower limit of mean diameter is 3 pm, unless A v is known by other means. 

(3) The spreads among the distributions are acceptably small. The spreads 
in sm are comparable to those in sp; the spreads in angle 9 assuming Mie 
diffraction are smaller than the spreads assuming Fraunhofer diffraction. 


Multiple Scatter Correction 

A model of multiple scatter was derived by reference 18, based on Fraun- 
hofer diffraction. Correction factors were presented for a Rosin-Rammler dis- 
tribution that determine diameter D and the dispersion parameter 6 from 
the measured D m , S m , and transmission factor x. The correction factors are 
Cp = D/D m and C$ = S/6 m . An exact correction like that obtained for Mie 
scatter would include seven distributions and four average diameters. In order 
to apply the method of reference 18 to this paper, two variables in the correc- 
tion equations of reference 18 were replaced by the variables used in this 
paper. 

(1) The dispersion parameter N or 6 was replaced by A v using the RR 
function in figure 9. 

(2) The obscuration OB = 1 - x was replaced by t. 

The two resulting equations and graphs of the equations of reference 18 
are shown in figure 16. 

The equations may be applied after the correction for Mie diffraction from 
the preceding section has yielded the corrected values of mean diameter and 
dispersion. The same correction is applied to all four diameters. 
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One limitation of the multiple scatter correction is the assumption of a 
transmission factor that is constant for every ray in the laser beam. This 
condition is not met with an axisymmetric spray where the transmission faactor 
is unity at the edge of the spray envelope, and much smaller near the center- 
line. To minize this source of error the laser beam cross-section should be 
adjusted to lie within the spray envelope. An alternative technique is first 
to use a fixed large diameter laser beam tc determine the mean diameter then to 
make the beam narrow in order to determine the transmission factor that can 
be used to determine the multiple scatter correction. 


ERROR ANALYSIS 

Gaussian Approximation cf Scattered Light 

In figure 17 the light scattered by perticles and transmitted by the scan- 
ning slit has a maximum weighted intensity Iw,max an 9 'e ©i , and a 
weighted intensity I w 2 at a fixed fraction of Iw.max a t angle 02- 
Weighted intensity I w ’ will be taken as irtensity I multiplied by 9. 

Using the Gaussian approximation of scattering (fig. 4) for the intensity 
I, the multiplied intensity is 

I w = 91 = 9 exp:-(aD9) 2 ] (14) 

where a = 0.57 ir/X , aD/0.57 = a, and 9 ; s in radians. 

The intensity ratio Iw,2^w,l a t ancles 02 and 9i is 

I w. 2 /I -.n,ax- < V 9 l ,e ’ t ' [-<*t» 2 ( 9 2 - «?)] n5) 

Angle Q] can be eliminated from equation (15) because ©i is the angle at 
which I w is a maximum. At this maximum 

(16) 


(17) 
is used as 

D varies 

inversely with © 2 ; however, an explicit solution for D is not possible. 

Differentiating equation (17) gives tie error in D caused by error in 
intensity measurement: 

dD _ (dI w,2 /I w,2 > ~ ( iI w,max /I w,max ) (18) 

D 1 - 2(a)0 2 ) 2 


0 2 = O.5(aD) 2 and I w max = <aD f^)“ 

Thus, equation (15) becomes 

Iw,2'Iw,max = aD0 2 exp [-(aD0 2 ) 2 

Figure 4 may be used to apply to the scann ng-slit method if (a02) 2 
the abscissa and (Iw.2^w,max ) ^ aD9 2 V 2 ^ s usec * as the ordinate. 

ion (17) shows that, at any intensity ratio, the diameter 
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when 02 is measured with neglible error. Because the various sources of 
extraneous light are independent of I , we assume that the random error in 
intensity measurement is constant and not dependent on I, and that the var- 
iances (dl w>max ) 2 and (dl w> 2 > 2 may each be set equal to (dl w ) 2 . 

Then equation (18) can be written as 


dD/D 

dI w /! 


w 


1 I <I w.max^ 1 




1/2 


1 - 2(aD0 2 )‘ 


(19) 


Equations (17) and (19) may be treated as parametric equations, with the 
common parameter aP 02 , so that the error ratio on the left side of equa- 
tion (19) may be represented as dependent solely on the ratio Iw,2^w,max- 
For 02 > 0i , the error ratio plotted in figure 18 shows a broad minimum 
between intensity ratios of 0.15 and 0.80. For 02 < 0] the error ratio plot- 
ted in figure 18 shows a broad minimum between intensity ratios of 0.15 and 
0.80. For 02 <0i, there is a narrower minimum between intensity ratios of 
0.55 and 0.85. The value of 0.85 is chosen to measure DvO.75 because it gives 
the largest measurement angle and thereby minimizes the effects of extraneous 
light caused by diffraction by the test section aperture. 

The minimum error is comparable to the minimum error of the two methods 
that were compared in reference 3, but with the important advantage that parti- 
cle diameter does not appear in figure 18. 


EXPERIMENTAL SIMULATIONS 

Confirmation of theoretical calculations appears possible by optical- 
laboratory experiments in which the spray and the background are simulated by 
glass plates carrying a dispersion of microspheres. In the experiments to be 
described, the measurements were made with a commercially available digital 
instrument designed for electrical waveform acquisition, storage, computation, 
and display, with output to a printer and a plotter. At the time of this test, 
a scanning pinhole was used instead of a scanning slit. 


Microsphere Calibration Plates and Mounting 

Polystyrene latex microspheres suspended in water are commercially avail- 
able, reference 23. A calibration plate is made by depositing the particles on 
a microscope cover slide preferably coated for minimum reflectance at the laser 
wavelength.. A typical coating is shown in figure 19 for 100 particles on 
both sides of the slide. The transmission factor of the particles on the slide 
is 0.85, or 0.92 for each side. To obtain a lower transmission factor, two 
slides are mounted 1 cm apart in a protective housing. To divert light that is 
multiply reflected by the slide surfaces, away from the detector, the two 
slides are inclined at an angle of 0.01 rad to each other, and the assembly is 
inclined at 0.01 rad or more to the plane-parallel test section windows. The 
two test section windows are also inclined at 0.01 rad to each other. Each 
slide is selected to have sides parallel to each other to 0.0001 radian as 
determined by inspection of interference fringes with a sodium vapor lamp. 
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A spray chamber to deposit particles cn a microscope slide is shown in 
figure 20. An artists airbrush spray-gun atomizes a suspension of microspheres 
in a liquid that has the same density as the particles. The liquid mixture is 
approximately 10 cc propanol, 6 cc Freon cleaning fluid, and 1 cc water con- 
taining the particles. Enough particles are included to produce a clouded mix- 
ture when viewed against a white background. The air brush fluid flow rate is 
adjusted by observing the spray-scattered light from the two microscope illu- 
minators. The upper illuminator shows the existence of spray. The lower illu- 
minator shows the degree of evaporation of the spray. The spray falling on 
the microscope slide must be sufficently evaporated to avoid formation of a 
liquid film on the slide, which will cause the particles to migrate and form 
clumps. Insufficient flow of spray will tend to produce dry particles or dust 
that does not accumulate on the slide. The concentration of particles on the 
slide is limited to that shown in figure IS. An attempt for greater particle 
concentration causes particles to accumulate on top of those already deposited. 


Analytical Correction f o • Extraneous Light 

In figure 21(a), a collimated laser beam with intensity I 0 is incident 
on a test section containing no spray. Sone of the light transmitted through 
the test section is not scattered. It has an intensity I 0 tq where x 0 is the 
transmission factor of the test section, including test section windows. 

A fraction of the light is forward-scattered with an intensity I 0 x 0 k o 
that is a fraction k 0 of the intensity of the axially-transmitted beam. The 
scattered-light intensity is usually less than one percent of the axial 
i ntensi ty . 

Separate detectors are used for the axial beam and for the scattered beam; 
these have calibration factors Cb and c e , respectively, each representing 
intensity per unit detector voltage. Equations (20) and (21) give the two 
detector voltages: E 0 for the axial beam and e o (0) for the beam scattered at 
angle 9. 


c a e 0 (9) = Iot 0 k 0 (9) (20) 

c b^o = Ic T o ^1 ) 

In figure 21(b) spray has been added lo the test section. Both the scat- 
tered and unscattered-1 i ght intensities given in equations (20) and 21) are 
equally attenuated by the spray particles with a transmission factor x s ; the 
complementary part 1 - x s is scattered by the spray in all directions. In 
the forward direction, the presence of spr<>y adds to the scattered light an 
intensity I o x o k s (0). The detector voltages are now e m (0) and E m (0) given 
by equations (22) and (23) as 

Ca6m (e) = ^o T o" c s^o + Io T o^s^®^ (22) 

c b^m = Io 1 o T s ^ 3 ) 

The desired quantity is l o xoks(0) for the spray alone. This is obtained by 
combining equations (20) through (23) to obtain 
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= c a e m (e) - x s c a e o (e) 


(24) 


where the spray transmission factor is 

t s = E m /E 0 (25) 

The desired intensity I o x o k s (0) alone would produce In the absence of back- 
ground, a voltage e s (9) is given by: 


l o x o k s (0) = c a e s (0) (26) 

Substitution of equation (26) in (24) gives the voltage distribution for the 
spray, corrected for background light, 

e s (0) = e m (0) - x s e o (0) (27) 

This is the general equation for correction of background light by subtraction 
of functions of 0. 

The validity of equation (27) was determined by a test with the micro- 
sphere calibration plates. If two plates are put in the test section 
(fig. 21(c)), one plate can be treated as an added cause of background light, 
and the other plate treated as spray. The intensity due to background can thus 
be made greater than the intensity due to spray in order to to exaggerate the 
needed correction. Both plates were first measured separately. For the back- 
ground plate (plate 1, subscript 1), the detector voltages analogous to equa- 
tions (22) and (23) with subscript 1 in place of subscripts s and m, are 
given by: 


c a ei (0) — I qx qx ^ k^ ( 0 ) + lQXQk](0) (28) 

c b E l = ^o^l <29) 

For the second plate (plate 2, subscript 2) alone, equations (28) and (29) 
change to 

c a e2(9) = Io^o T 2^o^®^ + Io T o^2^®^ (28b) 

c b E 2 = I o T o T 2 (29b) 

With both plates in the test section (fig. 21(c)) the intensities shown in 
the figure produce voltages given by: 

c a e 3(9) = Io T o x l T 2^o^®^ + ^o T o x 2 T 2^1 <6) + Io x o x 1^2^®^ (30) 

and 


c b E 3 = I o x o x l x 2 <31 > 

The desired intensity is represented by I o x o k2(0). Analogous to equa- 
tion (26), a voltage that would correspond to this intensity alone, if there 
were zero background, is represented by e s (0) where 

c a e s (0) = I 0 x 0 k 2 (e) (32) 
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Combination of equations (28), (30), and (12) yields 

e s (0) = Ce 3 (©) - t 2 ei(0)]/Ti (33) 

where 

tl = El/E 0 '2 = E 2 /E 0 (34) 

A check on the experiment is that the product xt . t 2 should be equal to 

e 3 /e 0 - 


Experimental Results 

Measurements were made with calibration plates bearing 25- and 50-pm 
microspheres to simulate the background an! the spray. The measured transmis- 
sion factors of these plates were 0.650 an 1 0.661 respectively. Figure 22(a) 
shows the voltage distributions e o (0) with no plates, e i (©) with the 25-pm 
plate alone, e 2 (0) with the 50-pm plate alone and, e3<©) with both plates. 

The voltage ratio E3/E0 = 0.423 with both plates is in good agreement with 
the product q. t 2 = 0.650 x 0.661 = 0.419. 

Figure 22(b) shows the true distribution e s (0) for the 25-pm plate as 
determined by equation (33), when the 50 pm plate was used to simulate the 
background. A comparison curve for the 25-pm plate alone, using equation (27) 
shows a smoother distribution and presumably is more accurate. There is good 
agreement between the curves; the background is thus only one ninth of the 
total signal. 

A second test of the equations is a comparison of the measurements of 
e 2 (0) using plate 2 alone with the calculation of e 2 (0) from measurements 
made where both plates are used. For the latter condition, equations (20), 
(28), and (30) give 


qe 2 (0) = xix 2 e o (0) - x 2 ei (0) + 63(0) (35) 

The two curves for e 2 (0), shown in figure 22(c), show good agreement. 

The correction for extraneous light i; in error in the presence of tur- 
bulence or other causes of refraction induced by the gas flow in the test 
section. The introduction of spray may change the the magnitude of these 
refraction effects, as shown by reference 5. For example, when our tests 
showed small-angle refraction by turbulence, the addition of spray reduced the 
refraction, so that the correction for extraneous light at small refraction 
angles was excessive. 


CONCLUSION 

A scanning optical system and a technique of data analysis have been 
described that permit determination of average particle diameters in a spray 
of droplets over a wide range of diameter distributions and of diameter disper- 
sions when there is a large aggregate of droplets. Mie computation shows a 
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diameter less than 3 pm cannot be accurately measured. Measurement of the 
angle at which the scattered-light intensity is 73 percent of the maximum 
intensity suffices to yield the volume median diameter. Measurement of the 
angle at which the scattered-light intensity is 42 percent of the maximum 
intensity suffices to yield the Sauter mean diameter. However, measurement of 
both of these angles and of the corresponding intensities at these angles also 
yields the dispersion of the volume size distribution. Measurement of the 
angle at which the scattered-light intensity is 85 percent of the maximum 
intensity yields the diameter at 0.75 of the cumulative volume. This diameter 
may be useful to measure the relative amount of large diameter particles. 
Measurement of the angle at which the scattered light intensity is 24 percent 
of the maximum intensity yields the volume/diameter mean D 31 . This diameter 
may be useful when compared with D 32 to better characterize the spray. The 
effect of extraneous light such as that introduced by the test section aperture 
may be reduced and corrected. The effect of multiple scatter may be corrected. 
Measurement of the spray transmission factor is needed to make both correc- 
tions. The effect of Mie scatter may be corrected. Use of a light source, 
like a laser, that produces a Gaussian distribution of off-axis irradiance 
improves accuracy. 

A method of using latex microspheres deposited on glass plates to simulate 
the spray and the background has been described. It provides a means of cali- 
brating the instrument and testing the measurement technique. Cross checks 
have indicated that the simulation is reliable. 
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APPENDIX - SYMBOLS 


A 

B 

Ca 

Cb 

C 

D 

D 

D 0 

°20 

d 30 

d 31 

d 32 

Dm 

D v0. 5 
DvO. 75 
D v0.1 
d v0.9 

Db 

Df 

D m,s 

Dp 

D u 

Dfi. 

d 

E 


0.57 

constant in distribution fun : t i on 

Constant in distribution function 

calibration constant for off -axi s-1 ight detector 

calibration constant for axiil-light detector 

constant in distribution function 

diameter of particle on aperture 

diameter at which volume distribution function is a maximum 

reference diameter 

diameter at mean of surface- area 

diameter at mean of volume 

diameter at mean of ( vol ume/ di ameter) 

Sauter mean diameter (at mean of volume/area) 
measured diameter 

volume median diameter, cumulative volume fraction 0.5 

diameter at cumulative volume fraction 0.75 

diameter at first decile of droplet volume distribution 

diameter at ninth decile of droplet volume distribution 

test section diameter 

diameter with Fraunhofer diffraction 

diameter with microsphere calibration 

particle diameter 

upper limit of particle diameter 

lower limit of particle dianeter 

aperture diameter 

voltage indicated by axial-light detector (a) 



e 


e 0- E 0 
e m> E m 
e s 

Eo to E 3 
eo to e 3 


f 

G 

H 


HO 

I 


io 

E max 

tw 

tw.max 
tw.max , tw2 
Jl 


K 

k 

k s 

k 0 

k] ,k 2 
N 


n 

P 

SF 

S M 


voltage indicated by scattered-light detector (b) 

value of e,E when there is no spray 

value of e,E where there is spray 

value of e for spray alone 

value of E under conditions 0 to 3 

value of e under conditions 0 to 3 

focal length of objective lens 

gain in intensity 

i rradi ance 

irradiance in axial direction 

radiant intensity 

radiant intensity of axial beam 

maximum radiant intensity 

weighted intensity, 19 

maximum weighted intensity 

weighted intensity at angle 9 ] ,02 

Bessel function of first kind 

scattering coefficient for single particle 

scattering coefficient 

scattering coefficient for spray 

scattering coefficient for empty test section 

scattering coefficient for plate 1, plate 2 

number of particles 

number fraction of particles 

radiant power 

slope parameter by assuming Fraunhofer diffraction 
slope parameter by assuming Mie diffraction 
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V 


Vo 

x 

*u 

y 

z 

a.9 

A V 

Ay ,m 
8 

*>m 

SA 

e 

e 

01 

02 
0F 
0M 

®m,s 

X 

Subscripts 

0 

s 

1 

2 

3 


volume of particles 
total volume of particles 
D/D 
D u /D 

ray displacement at focal distance f (fig. 1) 

dV/dD)/(dV/dD)p 

beam spread parameter 

dispersion of volume distribution 

measured A v 

parameter in distribution function 
measured 8 

area of scanning aperture (figs. 1 and 2) 
base of Napierian logarithms 
angular direction of forward -scattered ray 
angle at which maximum weighted Intensity occurs 
angle at which weighted intensity is measured 
angle assuming Fraunhofer diffraction 
angle assuming Mie diffraction 
angle with microsphere calibration 
wavelength of radiation 

transmission factor of test section 
transmission factor of spray 
transmission factor for plate 1 
transmission factor for plate 2 
transmission factor for plates 1 and 2 
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TABLE I. - NUMBER OF PARTICLES AT CONSTANT 
TOTAL VOLUME, AND A UNIT ORMLY-ILLUMINATED 
APERTURE DIAMETER OF 10 4 |juh WITH RESULT- 
ING TRANSMISSION FACTOR t AND 
INTENSITY GAIN G 


Diameter, 

pin 

N 

~ 

G 

IP - 

1 

0. C( 80 

1.01 ^ 

10* 

1(H 

.f-19 

1 .11x10;* 

10 1 

10 6 

.135 

3.20xl0 6 


TABLE II. - COORDINATES OF POINTS OF SMALLEST DEPENDENCE 


Mean 

diameter 

Scanni ng SI i t 

Scanning Pinhole 

, max 

a9 

A ( a0 ) / { a8 ) 

i/W 

«e 

A(a0) / {a9) 

D wfi c 

0.73 

2.10 

iO .08 

0.13 

2.64 

£0.05 

_vu . j 

Doo 

.42 

2.56 

£ . 16 

.04 

3.28 

£.12 


.85 

.91 

£.05 

.50 

1.75 

£.12 

_vu . t -J 
°31 

.24 

2.93 

£.25 

no acceptable point 


25 













FIGURE 6. - ENVELOPE OF RELATIVE- INTENSITY PEAKS OF FRAUNHOFER DIFFRAC- 
TION. WAVELENGTH OF RADIATION, X, 1 MM: PARTICLE DIAMETER, D; 

d 0 * 10^ UK- 
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RAY DEVIATION ANGLE, 9, DEG 


FIGURE 7. - ENVELOPE OF RELATIVE WEIGHTED- INTENSITY PEAKS OF FRAUNHOFER 
DIFFRACTION. WAVELENGTH OF RADIATION, X, 1 PARTICLE DIAMETER, D; 
d 0 - 10* 1 um. 



DROP iil/WETER. D. (IN 


FIGURE 8, - VOLUME DISTRIBUTION OF SIX FUNCTIONS WITH PEAK AMPLITUDE 
If = 50 pM, DISPERSION A v = 1.35. 
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RELATIVE INTENSITY, I/I^ OR VI W#MAX 



(3) WITH VOLUME MEDIAN Dy 0 5 . 




Jt. VU./3 

(b) WITH SAUTER MEAN D 32 . Cc) WITH 75 . 

FIGURE 12. - ENVELOPES OF SCATTERED- LIGHT DISTRIBUTION FOR SIX DISTRIBUTION FUNCTIONS AND FOR 
THREE VALUES OF VOLUME DISPERSION. 
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RELATIVE INTENSITY, I/I 










DISPERSION 



SLOPE PARAMETER, $p 


FIGURE 13. - CORRELATION BETWEEN VOLUME DISPERSION AND SLOPE PARAMETER FOR ALL 
DISTRIBUTION FUNCTIONS. 


APPROXIMATION FORMULA: 
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HOWEVER, IF (DyQ 5 <4.0 AND s^Sp >1.9), THEN S^/Sp = 1.9 
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FRAUNHOFER 
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DIAMETER, D^^, UM 

FIGURE 14. - SLOPE OF BEAM SPREAD CURVE FOR MIE SCATTER WITH WATER. 
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PARAMETER RATIO (a0) M /(a0) F 






PARAMETER RATIO <aO) M /(ae>F 



DIAMETER/ D 31 , pH 
(d) DIAMETER, D 3V 
FIGURE 15. - CONCLUDED. 





APPROXIMATION FORMULA: 


C D = 1.025 - (A v - 0.57)0.072 + 0.076 In T) 
WHERE C D > 1.0 AND T > 0.025 



(a) DIAMETER CORRECTION Cp. 



TRANSMISSION FACTOR. T 
(b) DISPERSION CORRECTION C A . 

FIGURE 16. - CORRECTION FACTOR FOR MULTIPLE SCATTER. 



0 ©1 0 2 

RAY DEVIATION ANGLE. 0 

FIGURE 17. - TWO-ANGLE ME \SUREMENT WITH SCANNING SLIT. 
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ERROR RATIO, (dD/D)/«J I/I, 



INTENSITY RATIO, V/Vl 
FIGURE 18. - ERROR RATIO, EQUATION (19). 



FIGURE 19. - MICROSPHERE CALIBRATION PLATE. 



f^-AIR IN 



WoW™ + I o T o lc s (0> 
Vo T s 

TEST 

SECTION 

(b) WITH SPRAY. 




FIGURE 20. - SPRAY CHAMBER FOR MICROSPHERE DEPOSIT ON MICROSCOPE FIGURE 21. - INTENSITY COMPONENTS TRANSMITTED BY TEST SECTION. 

SLIDE. 
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